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Quantum many-body systems divide into a variety of phases with very different physical prop- 
erties. The question of what kind of phases exist and how to identify them seems hard especially 
for strongly interacting systems. Here we make an attempt to answer this question for gapped in- 
teracting quantum spin systems whose ground states are short-range correlated. Based on the local 
unitary equivalence relation between short-range correlated states in the same phase, we classify 
possible quantum phases for ID matrix product states, which represent well the class of ID gapped 
ground states. We find that in the absence of any symmetry all states are equivalent to trivial 
product states, which means that there is no topological order in ID. However, if certain symmetry 
is required, many phases exist with different symmetry protected topological orders. The symmetric 
local unitary equivalence relation also allows us to obtain some simple results for quantum phases 
in higher dimensions when some symmetries are present. 

PACS numbers: 



I. INTRODUCTION 

For a long time, we believed that Landau symmetry- 
breaking theory 1 ' 2 describes all possible orders in materi- 
als, and all possible (continuous) phase transitions. How- 
ever, in last twenty years, it has become more and more 
clear that Landau symmetry-breaking theory does not 
describe all possible orders. For example, different frac- 
tional quantum Hall (FQH) states 3,4 all have the same 
symmetry. Thus it is impossible to use symmetry break- 
ing to characterize different FQH states. 

If Landau symmetry breaking theory is not enough, 
then what should we use to describe those new states 
of matter? It turns out that we need to develop a to- 
tally new theory, to described the new types of orders - 
topological/quantum order 5 ' 6 - that appear in the FQH 
states and the spin liquid states. 

In Ref. 7, a systematic understanding of a large class of 
topological orders in strongly correlated bosonic systems 
without symmetry has been developed based on string- 
net condensations. In Ref. 8, the string-net classification 
of topological orders was generalized, based on local uni- 
tary (LU) transformations. In Ref. 6,9,10, topological 
orders with symmetry are studied using projective sym- 
metry group and tensor network renormalization. But 
so far, we still do not have a complete classification of 
topological orders for interacting systems. 

Recently, for non-interacting gapped fermion systems 
with certain symmetries, a complete classification of 
topological phases has been developed based on the K- 
theory. 11,12 A generalization of the free fermion result to 
interacting cases has been obtained for ID systems. 13 

For ID bosonic systems, Ref. 14 studied quantized 
Berry phases for spin rotation symmetric systems. 
Ref. 15 studied the entanglement spectrum and the sym- 
metry properties of matrix product states. Using those 
tools, they obtained some interesting results which are 
special cases of the situation considered here. 



In this paper, we will apply the approach used in 
Ref. 8,17 to ID strongly correlated systems. We will com- 
bine the local unitary transformation with the symmetry 
properties of matrix product states, and try to obtain a 
complete classification of all gapped phases in ID quan- 
tum spin systems with certain simple symmetries. We 
find that 

(a) if there is no symmetry ( including translation sym- 
metry), all gapped ID spin systems belong to the 
same phase. 

(b) for ID non-translation invariant (NTI) spin sys- 
tems with ONLY an on-site symmetry described by 
a group G, all the phases of gapped systems that 
do not break the symmetry are classified by the 
equivalence classes in the second cohomology group 
H 2 (G, C) of the group G, provided that the physical 
states on each site form a linear representation of 
the group G. 

(Note that the equivalence classes in H 2 (G,C) classify 
the types of projective representations of G over the field 
C of complex numbers. Appendix C gives a brief intro- 
duction on projective representation and the second co- 
homology group.) In certain cases where G has infinitely 
many ID representations, for example when G = U(l), 
further classification according to different ID represen- 
tations exist. 

But quantum states are defined only up to global 
change of phases, therefore the symmetry operations of 
group G only needs to be represented by operators u{g) 
which satisfy 

u{gx)u{g 2 ) = e^^u( gi g 2 ) (1) 

for any group element <?i and g 2 . Such operators form a 
projective representation of group G. When we consider 
this general case, we find that the classification result 
remains the same as with linear representations, that is, 
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(c) for ID non-translation invariant (NTI) spin sys- 
tems with ONLY an on-site symmetry described by 
a group G, all the phases of gapped systems that 
do not break the symmetry are classified by the 
equivalence classes in the second cohomology group 
H 2 (G, C) of the group G, provided that the physical 
states on each site form a projective representation 
of the group G. 

In certain cases where G has infinitely many ID repre- 
sentations, further classification according to different ID 
representations exist. 

Applying these general results to specific cases allows 
us to reach the following conclusions: First, result (a) 
means that there is no non-trivial topological order in ID 
systems without any symmetry; Using result (b), we find 
that NTI integer-spin chain with only on-site 5*0(3) spin 
rotation symmetry can have two and only two different 
phases that do not break the SO (3) symmetry. Result 
(c) implies that NTI half-integer-spin chain with only on- 
site SO(3) spin rotation symmetry (which is represented 
projectively) also have two and only two gapped phase 
that does not break the SO (3) symmetry; We note the 
cyclic Z„ group has no non-trivial projective representa- 
tions, thus NTI spin chain with only on-site Z ra symmetry 
can have one and only one gapped phase that does not 
break the Z„ symmetry; The U(l) symmetry group has 
no non-trivial projective representation either, however, 
due to the special structure of the group of ID represen- 
tations of U(l), NTI spin chain with only on-site U(l) 
symmetry can have three and only three gapped phase 
that does not break the U(l) symmetry. 

We also considered systems with translation invari- 
ance (TI) and correspondingly many results have been 
obtained. 

(a) if there is no other symmetry, all gapped TI systems 
belong to the same phase. (This has been discussed 
as the generic case in Ref. 17). 

(b) for ID spin systems with ONLY translation sym- 
metry and an on-site symmetry described by a 
group G, all the phases of gapped systems that do 
not break the two symmetries are labeled by the 
equivalence classes in the second cohomology group 
H 2 (G,C) of the group G and different ID repre- 
sentations a{G) of G , provided that the physical 
states on each site form a linear representation of 
the group G. 

Similar to the NTI case, we should consider projective 
representations of G at each site. However, we find that, 

(c) there is no translation invariant gapped ground 
state symmetric under on-site symmetry of group G 
which is represented projectively on the state space 
at each site. 

In particular, we can show that the SO(3) spin ro- 
tation symmetric integer spin chain has two different 



gapped TI phases: 18-20 the spin-0 trivial phase and the 
Haldane phase 21 , while translation invariant 5*0(3) sym- 
metric half-integer spin chain must either be gapless or 
have degeneracy in ground space due to broken discrete 
symmetries 22,23 . On the other hand, the SU(2) sym- 
metric spin chains where on-site degrees of freedom con- 
tain both integer-spin and half-integer-spin representa- 
tions have only one gapped TI phase. We also show that 
the spin chain with only translation and parity symmetry 
(defined as exchange of sites together with an on-site Z-i 
operation) has four different gapped TI phases. 18-20 

For systems with time reversal symmetry, we find 
that NTI time reversal symmetric systems belong to two 
phases while TI time reversal symmetric phases on in- 
teger spin systems have two phases and those on half- 
integer spin systems are either gapless or have degeneracy 
in ground space. 

The paper is organized as follows: section II gives a de- 
tailed definition of gapped quantum phases and explains 
how that gives rise to an equivalence relation between 
gapped ground states within the same phase; section III 
shows that short-range correlated matrix product states 
represents faithfully ID gapped ground states and hence 
will be our object of study; section IV discusses the situ- 
ation where no symmetry is required and found no topo- 
logical order in ID; section V gives the classification of 
phases for ID systems with certain symmetries, for ex- 
ample on-site symmetry and time reversal symmetry. It 
classifies phases in translational invariant systems, and 
furthermore in systems where translational invariance is 
present together with other symmetries such as on-site 
symmetry, parity symmetry and time reversal symmetry; 
section VI generalizes some simple ID results to higher 
dimensions; finally in section VII we summarize our re- 
sults and conclude this paper. 



II. DEFINITION OF QUANTUM PHASES 

To obtain the above stated results, we need to first 
briefly discuss the definition of quantum phases. A more 
detailed discussion can be found in Ref. 8. Quantum 
phase describes an equivalence relation between quantum 
systems. Systems we consider live on a n-dimensional 
lattice and interactions are local (with finite range). A 
gapped quantum phase is usually defined as a class of 
gapped Hamiltonians which can smoothly deform into 
each other without closing gap and hence without any 
singularity in the local properties of ground state. Such 
an equivalence relation between Hamiltonians can be 
reinterpreted as an equivalence relation between ground 
states, as discussed in Ref. 8,17: Two gapped ground 
states belong to the same phase if and only if they are 
related by a local unitary (LU) transformation? 1 As LU 
transformations can change local entanglement structure 
but not the global one, states in the same phase have 
the same long range entanglement (LRE) and hence the 
same topological order. 5,24 States equivalent to product 
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states have only short range entanglement (SRE) and 
hence trivial topological order. All the states with short 
range entanglement belong to the same phase while states 
with long range entanglement can belong to the different 
phases. These considerations lead to the phase diagram 
as illustrated in Fig. 1(a), for the class of systems with- 
out any symmetry requirement. 

A LU transformation U can either take the form of 
finite time evolution with a local Hamiltonian 



(2) 



where T denotes time ordered integral and H(g) is a sum 
of local Hcrmitian terms. Or U can take the form of a 
constant depth quantum circuit 



(3) 



which is composed of R layers of unitaries and the 's 
within each layer k are local and commute with each 
other. 8 These two forms of LU transformations are equiv- 
alent to each other and we will mainly take the quantum 
circuit form for the discussion in this paper (the time evo- 
lution form is used for discussion of translation invariant 
systems) . More generally, we will consider equivalence re- 
lation between states defined on different Hilbert spaces 
and hence we allow a broader notion of unitarity (which is 
called generalized LU transformation in Ref. 8 and corre- 
sponds to the disentanglers and isometries in MERA 25 ). 

(k) 

Each LU operation U ik we consider in the quantum cir- 
cuit will act unitarily on the support space of the reduced 
density matrix of the region with (U^)^U^ = I of the 

original Hilbert space and U^ k \u^)^ = I of the sup- 
port space of reduced density matrix. While the total 
Hilbert space might change under such an operation, the 
entanglement structure of the state remains intact and 
hence the state remains in the same phase with the same 
topological order. 

If the class of systems under consideration have further 
symmetry constraints, two Hamiltonians are in the same 
phase if they can be connected by a smooth path that 
stays within this symmetric region. Correspondingly, the 
equivalence relation between gapped ground states of the 
same phase needs to be modified: 8 If the class of sys- 
tems have certain symmetry, two gapped ground states 
belong to the same phase if and only if they are related 
by a LU transformation which does not break the sym- 
metry. Such an restricted equivalence relation leads to 
a phase diagram with more structure, as shown in Fig. 
1(b). First not all short range entangled states belong to 
the same phase, short range entangled states with dif- 
ferent symmetry breaking will belong to different phases. 
Those symmetry breaking phases with short range entan- 
glement (SB-SRE) are well described by Landau's sym- 
metry breaking theory 1,26 . 

Can all phases with short range entanglement be de- 
scribed by symmetry breaking? Landau's symmetry 
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FIG. 1: (Color online) (a) The possible phases for class 
of Hamiltonians H(gi,g%) without any symmetry restric- 
tion, (b) The possible phases for class of Hamiltonians 
-ffsymm(<7i, 52) with some symmetries. The shaded regions in 
(a) and (b) represent the phases with short range entangle- 
ment. 



breaking theory suggest that states with the same sym- 
metry always belong to the same phase, which implies 
that all phases with short range entanglement are de- 
scribed by symmetry breaking. However, this result 
turns out to be not quite correct. States which do not 
break any symmetry can still belong to different phases 
as well. 6 We will refer to the order in symmetric short 
range entanglement (SY-SRE) states as symmetry pro- 
tected topological order 19 . For example, in the presence 
of parity symmetry, the Haldane phase and the S z — 
phase of spin-1 chain belong to two different phases even 
though both phases have short range entanglement and 
do not break the parity symmetry. Also, in the presence 
of time-reversal symmetry, the topological insulators and 
the band insulators belong to two different phases. Again 
both phases have short range entanglement. 

For systems with long range entanglement, the phase 
diagram similarly divides into symmetry breaking (SB- 
LRE) and symmetric (SY-LRE) phases. The charge 
4e superconducting states 27 and the symmetric Z2 
states 6, 28 are examples of the SB-LRE phases and the SY- 
LRE phases respectively. 



III. ID GAPPED SPIN SYSTEMS AND 
MATRIX PRODUCT STATES 

Having defined the universality classes of phases as the 
equivalence classes of states under (symmetric) local uni- 
tary transformations, we would then like to know which 
phases exist, or in other words, to classify all possible 
phases in strongly correlated systems. Some partial clas- 
sifications have been discussed for strongly correlated sys- 
tems through string-net states, 7 and for free fermion sys- 
tems with certain symmetries through K-theory. ' In 
this paper, we would like to consider ID gapped strongly 
correlated spin systems both with and without symme- 
try, and try to classify all such systems whose ground 
state does not break any symmetry. (In other words, 
the ground state has the same symmetry as the Hamil- 
tonian.) 
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Completely classifying strongly correlated spin systems 
seems to be a hard task as in general strongly interact- 
ing quantum many-body systems are very hard to solve. 
However, the recent insight about describing ID gapped 
ground states of spin systems with matrix product state 
formalism 29 ' 30 provides us with a handle to deal with this 
problem. A matrix product states (MPS) is expressed as 

\<j>) = ]T Tr(A^A^...AfJ)\hi 2 ...i N ) (4) 

ii,i 2 ,...,ixr 

where ik = l...d with d being the physical dimension of a 
spin at each site, A 1 - J 's are D x D matrices on site k with 
D being the inner dimension of the MPS. It has been 
shown that matrix product states capture the essential 
features of ID gapped ground state, for example an en- 
tanglement area law 31 and finite correlation length 32 ' 33 , 
and provide an efficient description of such states 34 . On 
the other hand, generic matrix product states satisfying a 
condition called 'injectivity' are all gapped ground states 
of local Hamiltonians 29,30 . Therefore, studying this class 
of MPS will enable us to give a full classification of ID 
gapped spin systems. 

Now the question of what phases exist in ID gapped 
spin systems can be restated as what equivalence classes 
of matrix product states exist under LU transformations. 
Ref. 17 gave a specific way to apply such LU transforma- 
tions, which realizes a renormalization group transforma- 
tion on MPS that removes local entanglement and takes 
the states to a simple fixed point form. A partial classi- 
fication of MPS is also given in Ref. 17. In the following 
we will use this procedure to classify gapped phases of 
ID spin system, in particular the ID systems with vari- 
ous symmetries. We see that the possible phases in ID 
strongly correlated systems depend on the symmetry of 
the class of systems. 

First we will briefly review how the renormalization 
group transformation 1 ' is done. For the identification 
and optimal removal of local entanglement, a particularly 
useful mathematical construction is the double tensor 

eJU = x ( A S*)* of the MPS - Elk] unk i uel y 

determines the state up to a local change of basis on each 
site 30,35 , that is, if 

F W -V4 [fc] x (A [k] r-VB W x (R [k] Y 

i i 

(5) 

then A^ afj and B^^ are related by a unitary transfor- 
mation 

3 

Therefore, states described by A^ and B^ 1 have exactly 
the same entanglement structure, which is faithfully cap- 
tured in E^'l . A proof of this fact can be found in Ref. 35. 
For clarity, we present the proof in appendix A following 
the notation of this paper. 



Local unitary operations on MPS can be applied 
through manipulation of . Treat as a D 2 x D 2 ma- 
trix with row index aj and column index f3\. To apply 
a unitary operation on n consecutive sites, we combine 
the double tensor of the n sites together into 

£ = E [1] E [2] ...E [n] (7) 
and then decompose E into a set of matrices Aj's 

i 

Note Ai a p is determined up to a unitary transformation 

on i. The index i of A^ a g, up to an unitary transforma- 
tion, can be viewed as the combination of ii, i%, i n , the 
indices oi A [ ^ afs , A [ £ aj3 , ... A l ^ a/3 . Going from original 

indices i\, %2 1 i n to the effective index i corresponds 
to applying a unitary operation on the n-block and A^ 
describes the new state after operation. 

The unitary operation can be chosen so that local en- 
tanglement is maximally removed. E contains all the in- 
formation about the entanglement of the block with the 
rest of the system but not any detail of entanglement 
structure within the block. Hence we can determine from 
E the optimal way of decomposition into A which corre- 
sponds to the unitary operation that maximally removes 
local entanglement while preserving the global structure. 
To do so, think of E aj ^ x as a matrix with row index 
a/3 and column index 7%. It is easy to see that with 
such a recombination, E is a positive matrix and can be 
diagonalized 

i 

where we have kept only the non-zero eigenvalues Aj and 
the corresponding eigenvectors V\ a p- A\s then given by 

\ aP = \Ai*W (10) 

which are the matrices representing the new state. In 
retaining only the non-zero eigenvalues, we have reduced 
the physical dimension within the block to only those 
necessary for describing the entanglement between this 
block and the rest of the system. Local entanglement 
within the block has been optimally removed. 

Each renormalization step in the renormalization pro- 
cedure hence works by grouping every n consecutive sites 
together and then applying the above transformation to 
map AW, A^- 2 \..., AM to A. So one renormalization step 
maps the original matrices (A\ k ')(°} on each site to renor- 

malized matrices (A^)^ on each block. Repeating this 
procedure for a finite number of times corresponds to 
applying a finite depth quantum circuit to the original 
state. If the matrices reaches a simple fixed point form 
(A^)(°°) (up to local unitaries), we can determine from 
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it the universal properties of the phase which the original 
state belongs to. 

Such a renormalization procedure hence provides a way 
to classify matrix product states under LU transforma- 
tions by studying the fixed point (A^ that a state 
flows to. Two states are within the same phase if and 
only their corresponding fixed points states can be trans- 
formed into each other by (symmetric) LU transforma- 
tions. In the following we will apply this method to study 
short-range correlated matrix product states which faith- 
fully represent the class of ID gapped ground states. 

The short-range correlation is an extra constraint on 
the set of matrix product states that we will consider. 
Not all matrix product states describe gapped ground 
states of ID spin systems. In particular, ID gapped 
ground states all have finite correlation length 33 for equal 
time correlators of any local operator, while matrix prod- 
uct states can be long-range correlated. The finite corre- 
lation length puts an extra constraint on MPS that the 
eigen-spectrum of E should have a non degenerate largest 
eigenvalue (set to be 1) (see appendix B). 29 ' 30 Therefore, 
we will assume this property of E in our following discus- 
sion and corresponding MPS will be called short-range 
correlated (SRC) MPS. 

This renormalization method is well suited for the 
study of systems without translational invariance, which 
we will discuss in detail. We will also make an attempt to 
study translational invariant systems with this method. 
While the full translational symmetry is reduced to block 
translational symmetry in the RG process, by studying 
the resulting equivalence classes for different values of 
block size n, we expect to obtain a more complete classifi- 
cation of translational invariant ID systems. Indeed, the 
classification result is further confirmed by using a trans- 
lational invariant LU transformation in the time evolu- 
tion form to study equivalence between TI systems. 

In the following sections, we will present our analysis 
and results for different cases. First we will consider the 
situation where no specific symmetry is required for the 
system. 

IV. NO TOPOLOGICAL ORDER IN ID 

When no symmetry is required for the class of system, 
we want to know what kind of long range entanglement 
exists and hence classify topological orders in ID gapped 
spin systems. We will show that 



All gapped ID spin systems belong to the same 
phase if there is no symmetry. 



In other words, there is no topological order in ID. This 
is similar to the generic case discussed in Ref. 17. 

To obtain such a result, we use the fact that gapped 
ID spin states 52 are described by short-range correlated 
(SRC) matrix product states. Then one can show that 
all SRC matrix product states can be mapped to product 



states with LU transformations and hence there is no 
topological order in ID. 

Consider a generic system without any symmetry (in- 
cluding translation symmetry) whose gapped ground 

state is described as an MPS with matrices Af^ that 
vary from site to site. Ref. 30 gives a 'canonical form' 

for the matrices so that the double tensor E 1 ' ^ , when 
treated as a matrix with row index and column index 

fix, has a left eigenvector Aa\ = Aa'<5 Q7 and correspond- 
ing right eigenvector Ajg* 1 ' = A^° +1 '^ x . Here A's are 

positive numbers and J2 a ^a = 1- ^a-y = 1 when a = 7 
and 8 ai = otherwise. 53 This eigenspace has the largest 
eigenvalue in E[ fe l 36 and is usually set to be 1. Note that 
the right eigenvector on site k is the same as the left 
eigenvector on site k + 1 and has norm 1, therefore when 
multiplying the double tensors together, this one dimen- 
sional eigenspace will always be of eigenvalue 1 . 

There could be other eigenvectors of eigenvalue 1 in 
EW. However, this will lead to an infinite correlation 
length 29,30 and hence not possible in ID gapped state. 
Therefore, for short-range correlated MPS, E^ must 
have a non-degenerate largest eigenvalue 1. When multi- 
plying the double tensors together, the remaining block 
of E^l will decay exponentially with the number of sites. 
This consideration is essential for determining the fixed 
point of the renormalization procedure when applied to 
the MPS, as shown below. 

Now we apply the renormalization procedure as dis- 
cussed in the previous section to remove local entan- 
glement from a general SRC MPS. Take block size n. 
The double tensor on the renormalized sites are given by 
(£[#])« = n fce .R-(E [fcl ) (0) , where k's are the n sites in 
block K. (E, again, is treated as a D 2 x D 2 matrix with 
row index a"f and column index fix-) 

After repeating the renormalization process a finite 
number of times, (E^)^ will be arbitrarily close to a 
fixed point form (E^l)^ 00 ) with non-degenerate eigenvalue 
1 and (EW)£3, X = AL fe | Af where AL^ 1 = \£ ] S ai and 
r[fe+i] _ x[fc+i] A 

Now we can decompose (E^^ 00 ) into matrices to find 
the fixed point state. One set of matrices giving rise to 
this double tensor is given by 

(41^ = ^^-^^^ (ii) 

i l ,i r = 1...D. Here we use a pair of indices (i l , i r ) to label 
the effective physical degrees of freedom on the renormal- 
ized site k, and (A^\ r )^°°^ is a set of matrices that defines 
the fixed-point MPS. It is clear from the form of the ma- 
trices that at fixed point every site is composed of two 
virtual spins of dimension D. Every virtual spin is in an 
entangled pair with another virtual spin on the neighbor- 
ing site \EPk : k+i) = J2iLi Af^^lM) an d the full many- 
body state is a product of these pairs. An illustration of 
this state is given in Fig. 2(upper layer). 
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FIG. 2: (Color online) Disentangling fixed point state (upper 
layer, product of entangled pairs) into direct product state 
(lower layer) with LU transformations. 



Obviously we can further disentangle these pairs by ap- 
plying one layer of local unitary transformations between 
every neighboring sites and map the state to a product 
state (Fig. 2, lower layer). 

Therefore, through these steps we have shown that all 
SRC matrix product states can be mapped to product 
states with LU transformations and hence there is no 
topological order in ID NTI system. 



SYMMETRY PROTECTED TOPOLOGICAL 
ORDER IN ID 



where 0(<7i,<72) = in a linear representation and 
#(<7i, S2) could take non-trivial value in a projective rep- 
resentation. A projective representation of a symme- 
try group is generally allowed in a quantum description 
of system because the factor e 16 ^ 1 ' 92 ' only changes the 
global phase of a quantum state but not any physically 
measurable quantity. Therefore, in our classification, we 
will consider not only the case of linear representation, 
but also in general projective representations. 

The on-site symmetry is the only symmetry required 
for the class of system. In particular, we do not require 
translational symmetry for the systems. However, for a 
simple definition of phase, we will assume certain uni- 
formness in the state, which we will define explicitly in 
the following. We will classify possible phases for differ- 
ent G when the ground state is invariant (up to a to- 
tal phase) under such on-site symmetry operations and 
is gapped (ie short-range correlated). Specifically, the 
ground state \4>l) on L sites satisfies 



u^{g)®...®u^{g)\^ L )=a L {g)\<j> L ) (13) 
where |ai(g)| = 1 are g and L dependent phase factors. 



If the class of systems under consideration has certain 
symmetry, the equivalence classes of states are defined 
in terms of LU transformations that do not break the 
symmetry. Therefore, when applying the renormalization 
procedure, we should carefully keep track of the symme- 
try and make sure that the resulting state has the same 
symmetry at each step. Due to such a constrain on lo- 
cal unitary equivalence, we will see that gapped ground 
states which do not break the symmetry of the system 
divide into different universality classes corresponding to 
different symmetry protected topological orders. We will 
first discuss the case of on-site symmetries in detail for 
non-translational invariant (NTI) systems, i.e., the sys- 
tem has only an on-site symmetry and no translation 
symmetry. Then we shall make an attempt to study 
translational invariant (TI) systems, with the possibility 
of having on-site symmetry or parity symmetry. Lastly, 
we shall consider the case of time reversal symmetry. 

A. On-site Symmetry 

A large class of systems are invariant under on-site 
symmetry transformations. For example, the Ising model 
is symmetric under the Z 2 spin flip transformation and 
the Heisenberg model is symmetric under SO (3) spin ro- 
tation transformations. In this section, we will consider 
the general case where the system is symmetric under 
v,(°\g) ® ... ® v,(°\g) with u^°\g) being a unitary rep- 
resentation of a symmetry group G on each site. The 
representation can be linear or projective. That is, for 
any g 1 ,g 2 G G, 

u{ 9l )u{g 2 ) = e^^u(g m ) (12) 



1. On-site Linear Symmetry 

First, let us consider the simpler case where uS°'(g) 
form a linear representation of G. at (g) is then a one- 
dimensional linear representation of G. Now we will try 
to classify these symmetric ground states using symmet- 
ric LU transformations and we find that 

Consider ID spin systems with ONLY an on- 
site symmetry G which is realized linearly, all 
the gapped phases that do not break the symmetry 
are classified by H 2 (G,C), the second cohomology 
group ofG, if H 2 (G,C) is finite andG has a finite 
number of ID representations. 

We will also discuss the case of 17(1) group which has an 
infinite number of ID representations. 

We will again assume that all gapped states can be rep- 
resented as short range correlated matrix product states. 
We will use the renormalization flow used before 17 to sim- 
plify the matrix product states and use the fixed-point 
matrix product states to characterize different equivalent 
classes of LU transformations, as two symmetric states 
belong to the same class if and only if their correspond- 
ing fixed-point states can be mapped to each other with 
symmetric LU transformations. 

In order to compare different equivalent classes under 
symmetric LU transformations, it is important to keep 
track of the symmetry while doing renormalization. 

First, in the renormalization procedure we group n 
sites together into a new site. The on-site symmetry 
transformation becomes tv-°\g) = (®u(°'(g)) n , which is 
again a linear representation of G. The next step in RG 
transformation applies a unitary transformation u?[ to 
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the support space of new site k. This is actually itself 
composed of two steps. First we project onto the sup- 
port space of the new site, which is the combination of 
n sites in the original chain. This is an allowed opera- 
tion compatible with symmetry G as the reduced den- 
sity matrix p n is invariant under u^(g), so the support 
space form a linear representation for G. The projection 
of u^(g) onto the support space P n u^°'(g)P n hence re- 
mains a linear representation of G. In the next step, we 
do some unitary transformation iy[ within this support 
space which relabels different states in the space. The 
symmetry property of the state should not change under 
this relabeling. In order to keep track of the symmetry of 
the state, the symmetry operation needs to be redefined 
as = wf ] P n u^(g)P n (wf ] y. After this re- 

definition, the symmetry operations (u^)^(g) on each 
new site k form a new linear representation of G. 

By redefining (W fe ')^(<?) at each step of RG transfor- 
mation, we keep track of the symmetry of the system. Fi- 
nally at the fixed point (ie at a large RG step i = R), we 
obtain a state described by (A^\ r ) ^ which is again given 
by the fixed point form eqn. (11). To describe a state that 
does not break the on-site symmetry, here (A^. r ) ^ is in- 
variant (up to a phase) under (v} k ')( R '(g) on each site k. 
Therefore 16 , 

Y, u Wr^ A % = *\!^9)N^{g)A%M [k] (g) 
j l j r 

N [k] (g) = M^ig) (14) 

must be satisfied with some invertible matrix (g) and 
M\k](g)- Here k labels the coarse grained sites and we 

{ R) 

have dropped the RG step label R (except in a) k -, (g)). 
Each coarse grained site is a combination of n R original 

( R) 

lattice sites and ofa (g) form a ID (linear) representation 
of G. 

Solving this equation we find the following results (see 
appendix D): 

(a) Nu^(g) and M^(g) are projective representations of 
G (see eqn. (D7)). Projective representations of G be- 
long to different classes which form the second cohomol- 
ogy group H 2 (G,C) of G. (For a brief introduction on 
projective representation, see appendix C). M^(g) and 
N[k\(g) corresponds to the same element ui in H 2 (G,C). 

(b) The linear symmetry operation (g) must be of the 
form (g)u^' l (g) (^u^' r (g) where u^' 1 and u^' r are 
projective representations of G and correspond to inverse 
elements ui and — cu in H 2 (G,C) respectively. a^(g) is 

a ID (linear) representation of G. u^' 1 and vJ- k ^ r act on 
the two virtual spins separately (see eqn. (D13)). 

Therefore, the fixed point state is formed by entangled 
pairs l-E-Pfc.fc+i) of virtual spins which are invariant, up 

to a phase (due to the non-trivial a, k ^ (<?)), under linear 

transformation u^' r (g) ® u^^'^g). 

Now we use the uniformness of the state and simplify 



our discussions. Specifically, we assume that aL '(g) docs 

( R) 

not depend on the site index k. Certainly, CKm (<?) does 
not depend on k if the state has the translation sym- 
metry. If the ID representations of G are discrete, then 
for weak randomness that slightly break the translation 
symmetry, (g) still does not depend on k. So we can 

drop the k index and consider aS R \g). 

Does different a^ R \g) label different symmetric 
phases? First, the answer is no if the number of ID rep- 
resentations of G is finite (as is the case for Z ni 50(3), 
etc) . Because we can always choose block size n properly 
so that a^ R \g) = 1 and the difference between symmet- 
ric states due to a^(g) disappears. In the case of U(l) 
group, there are infinitely many different ID represen- 
tations e %mB , labeled by integer m. For two states with 
positive mi, 777,2, we can always choose block size m2n R , 
m\n R respectively, so that the ID representations be- 
come the same. This is also true for negative mi, m 2 . 
But if mi, 77i2 take different sign (or one of them is 0), 
the ID representations will always be different no matter 
what blocking scheme we use. Therefore, U(l) symmet- 
ric states divide into three classes due to different ID 
representations that can be labeled by {+,0,—}. After 
these considerations, we will ignore the ID representa- 
tions oS R \g) in the following discussion. 

We find that the entangled pairs \EPk y k+i) of virtual 
spins in the fixed point state are exactly invariant un- 
der linear transformation u^' r (g) ® v} k+1 ^ L (g). The left 
virtual spin of each site forms a projective representa- 
tion of G corresponding to element ui in H 2 (G, C), while 
the right virtual spin corresponds to element — us. In 
appendix E, we will show that fixed point states with 
the same us can be related by a symmetric LU transfor- 
mation, while those with different us cannot. Therefore, 
the phases of SRC MPS that are invariant under lin- 
ear on-site symmetry of group G are classified by the 
second cohomology group H 2 (G,C). (When G — U(l), 
further division of classes due to different ID representa- 
tions of G exist. The equivalence classes are labeled by 
a e {+,0,-} and ui e H 2 (U(1), C).) 



2. On-site Projective Symmetry 

Due to the basic assumption of quantum mechanics 
that the global phase of a quantum state will not have 
any effect on the physical properties of the system, it is 
necessary to consider not only the linear representation 
of symmetry operations on the system, but also the pro- 
jective representations. For example, on a half-integer 
spin, rotation by 2tt is represented as — /, minus the iden- 
tity operator instead of I. Hence, the rotation symmetry 
50(3) is represented projectively on half integer spins. In 
order to cover situations like this, we discuss in this sec- 
tion systems with on-site projective symmetry of group 
G. 

Again, we consider the case when the ground state does 



not break the symmetry, i.e. u^°\g)<S> ... <S> {g)\4>L) — 
a L(g)\4>L) , where u^°\g) form a projective representation 
of group G corresponding to class u. Assuming uniform- 
ness of the state, we require that u) does not vary from 
site to site. 

But this can be reduced to the previous linear case. As 
long as H 2 (G,C) is finite and w has a finite order n, we 
can take block size n so that after blocking, the symmetry 
operation on the renormalized site u^(g) = (CsW-' (<?))" 
corresponds to nuj — luq in H 2 (G,C). Therefore, the 
state after one blocking step is symmetric under an on- 
site linear representation of group G and all the reasoning 
in the previous section applies. We find that the classi- 
fication with projective on-site symmetry is the same as 
linear on-site symmetry. That is 

Consider ID spin systems with ONLY an on- 
site symmetry G which is realized projectively, all 
the gapped phases that do not break the symmetry 
are classified by H 2 (G,C), the second cohomology 
group ofG, if H 2 (G,C) is finite andG has a finite 
number of ID representations. 

The U(l) group does not have a non-trivial projective 
representation and will not introduce any complication 
here. 

3. Examples 

Since G = Z„ has no non-trivial projective representa- 
tions, we find that 

All ID gapped systems with only on-site Z„ sym- 
metry belong to the same phase. 

For spin systems with only spin rotation symmetry, 
G = 50(3). 50(3) has two types of projective represen- 
tations described by H 2 (SO(3),C) = {0, 1}, correspond- 
ing to integer and half-integer spin representations. We 
find that 

For integer-spin systems, all ID gapped systems 
with only on-site 50(3) spin rotation symmetry 
have two different phases. 

Such a result has some relation to a well known result 37 
for NTI spin-1 Heisenberg chain 

i 

The model undergoes an impurity driven second order 
phase transition from the Haldane phase 21 to the random 
singlet phase 38 ' 39 as the randomness in Ji increases. 

For half-integer-spin systems, 50(3) is represented 
projectively on each site, yet the classification is the same 
as the integer case, we find that 

For half-integer-spin systems, all ID gapped states 
with only on-site 50(3) spin rotation symmetry 
have two different phases. 



Representative states of the two phases are nearest- 
neighbor dimer states, but with dimer between sites 2i 
and 2i + 1 in the first phase and between sites 2i — l and 
2i in the second phase. 

The projective representation of 50(3) on half-integer- 
spins form a linear representation of SU(2). If we think 
of the linear representation of 50(3) on integer-spins as a 
(unfaithful) linear representation of SU(2) and allow the 
mixture of integer and half-integer spins on one site, then 
the two phases of 50(3) merge into one. Therefore, sys- 
tems with only on-site SU (2) symmetry (which implies 
the mixture of integer and half-integer spins on each site) 
belong to one phase as we can map integer-spin singlets 
into half-integer-spin singlets without breaking the SU(2) 
symmetry (see appendix E). Such a procedure breaks 
down if 50(3) symmetry is required for each site as the 
direct sum of a linear representation (on integer-spin) 
and a projective representation (on half-integer-spin) is 
no long a projective representation for 50(3). 

In this way, we have obtained a full classification of 
the phases of gapped NTI ID spin systems with various 
on-site symmetry. 

B. Translation Invariance 

We have discussed the gapped phases of ID NTI sys- 
tems that have some on-site symmetries. In this section, 
we would like to discuss translation symmetric systems. 
We will consider those translation symmetric systems 
whose ground states are gapped and translation invariant 
(TI). 54 

1. {ni}-block TI LU transformations and b-phases 

To discuss the TI phases, we need to discuss the equiv- 
alence classes under TI LU transformations. However, it 
is hard to use quantum circuit to describe TI LU trans- 
formations. Thus, the quantum circuit formulation used 
in this paper is inconvenient to describe TI LU transfor- 
mations. However, in this section, we will first try to use 
the quantum circuit formulation of LU transformations 
to discuss the phases of TI gapped states. While we are 
able to identify many different phases in this way, we can 
not rigorously prove the equivalence of states within each 
phase using the non-translational invariant circuit. In or- 
der to establish this equivalence, we later apply the other 
formulation of local unitary transformation - a finite time 
evolution with a local Hamiltonian(Eqn.2) where TI can 
be preserved exactly and we confirm the classification 
result obtained with quantum circuits(see appendix G). 

Let us consider the LU transformations represented by 
quantum circuits which is formed by the unitary opera- 
tors on blocks of ni sites in the i th layer (see Fig. 3). We 
will call such LU transformations {nj}-block LU trans- 
formations. If the unitary operators on different blocks 
in the same layer of the quantum circuit are the same, 
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FIG. 3: A {ni}-block TI LU transformation described by a 
quantum circuit of three layers. Here n\ = 4, ni = 3, riz = 2. 
The unitary transformations Ui on different blocks in the i th 
layer are the same. 

we will call the LU transformations a {n.;}-block TI LU 
transformations. 

In this paper, we will try to use quantum circuit for- 
mulation of the LU transformations to discuss gapped 
translation symmetric phases. One way to do so is to 
use the equivalent classes of the {rii}-block TI LU trans- 
formations to classify the gapped translation symmetric 
phases. However, since the {ni}-block TI LU transforma- 
tions are different from the TI LU transformations, the 
equivalent classes of {ni}-block TI LU transformations 
are different from the equivalent classes of the TI LU 
transformations. But since the TI LU transformations 
are special cases of {ni}-block TI LU transformations, 
each equivalent class of {rii}-block TI LU transforma- 
tions are formed by one or more equivalent classes of TI 
LU transformations. 

Therefore, we can use the equivalent classes of {rii}- 
block TI LU transformations to describe the gapped TI 
phases, since different equivalent classes of {n^l-block TI 
LU transformations always represent different TI quan- 
tum phases. On the other hand, the equivalent classes of 
{rii }-block TI LU transformations may not separate all 
gapped TI phases. Sometimes, a single equivalent class 
of {ni}-block TI LU transformations may contain several 
different gapped TI phases. 

To increase the resolution of the {ni}-block TI LU 
transformations, we would like to introduce the block 
equivalent classes: two states belong to the same block 
equivalent classes, iff, for all values of rii, they can be 
mapped into each other through {ni}-block TI LU trans- 
formations. 

Clearly, each block equivalent class might still contain 
several different gapped TI phases. In this paper, we 
will call a block equivalent class a block phase, or b- 
phase. It is possible that the block equivalent classes 
are the same as the universal classes that represent the 
gapped TI phases. In this case the gapped TI b-phases 
are the same as the gapped TI phases, and we can study 
the gapped TI phases through block equivalent classes. 
In the following, we will first study the b-phases for ID 
strongly correlated systems with translation and some 
other symmetries and then confirm that gapped TI b- 
phases do coincide with gapped TI phases. 

To describe TI gapped states we will use TI MPS rep- 



resentation with site independent matrices. But how do 
we know a TI gapped state can always be represented by 
a MPS with site independent matrices? A non-uniform 
MPS in eqn. (4) can represent a TI state if the matrices 
satisfy, for example, 

A [k+i] = M -i A [k] M (16) 

for an invertible matrix M. 

It is proven in Ref. 30 that every TI state does have 
a TI MPS representation. Specifically, in the example 
considered above, we can transform the matrices so that 
they become site independent. Let us introduce 

By construction, and Af^ will represent the same 
MPS. We see that 

Alk+l] _ N -i A [k+i] N 

= N^M^AfMN^ (18) 
= N^M- 1 N [k _ 1] A [ ? ] N [k lMN [k+1] 

We see that if we choose 

N [k] = M-\ (19) 

we will have 

if +1 J=if]. (20) 

Hence reducing the original non-uniform representation 
to a uniform one. 

We also like to remark that by "TI gapped state" , we 
mean that there exists a TI gapped Hamiltonian on 
a lattice of L sites such that the TI gapped state is the 
ground state of Hl ■ We would like to stress that we do 
not need the above condition to be true for all values of L. 
We only require Hz, to be gapped for a sequence of lattice 
sizes {Li} that lim^oo Li = oo. In this paper, when 
we discuss a system of size L, we always assume that 
L belongs to such a sequence {Li}. This consideration 
is important if we want to include in our discussion, for 
example, boson systems with 1/2 particles per site which 
can only be defined for even system size L. 



2. Gapped TI b-phases coincide with gapped TI phases 

In the above discussion, we see that for a large block 
size, the related matrix [see (7)] E = eWE^I-EM is dom- 
inated by its largest eigenvalue. For translation invari- 
ant states, EW = E and E = E n . So the fixed-point E, 
and hence the gapped TI b-phase, is directly determined 
by the largest eigenvalue and the corresponding left and 
right eigenvectors of E. 
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Since a gapped TI b-phase is directly determined by E 
on each site, we do not need to do any blocking transfor- 
mation to understand the gapped TI b-phase. We can di- 
rectly extract its fixed-point tensor from E. This suggests 
that gapped TI b-phases coincide with gapped TI phases. 
Indeed, we can directly deform E to the fixed-point E by 
changing other non-largest eigenvalues to zero. Such a 
procedure allows us to deform the matrix Aj jQ| g of the 
initial MPS to the fixed-point matrix a a of the final 
MPS. Since the largest eigenvalue of E is not degenerate, 
the state remains short range correlated (and gapped) 
during the deformation. 29,30 (As we are writing up the 
above arguments, we learned that a similar method is 
used by Schuch et. al. 47 ) Also the state does not break 
the translation symmetry (unlike the n-block transfor- 
mations) and other symmetries during the deformation. 
This allows us to show that the gapped TI phases is char- 
acterized by the largest eigenvalue and the correspond- 
ing left and right eigenvectors of E. Thus gapped TI 
b-phases coincide with gapped TI phases. For details, 
see appendix G. 

In the following, we will use the {rti}-block TI LU 
transformations to discuss various gapped TI phases of 
ID systems. 

3. ID systems with only translation symmetry 



For ID systems with only translation symmetry, 
there is only one gapped TI phase. 

Such a result generalizes the earlier result that for ID sys- 
tems with no symmetry, there is only one gapped phase. 
To obtain the new result, we basically repeat what we 
did in section IV. The only difference is that the matri- 
ces representing the state now are site independent and 
in section IV we use {rii}-block LU transformations to 
reduce the ID NTI MPS, while to derive the new result, 
here we use {rij }-block TI LU transformations to reduce 
the ID TI MPS. 



4- ID systems with translation and on-site unitary 
symmetries 

Similarly, by repeating the discussions in section V A 1 
and section VA2 for {ni}-block TI LU transformations 
on the ID TI MPS, we can show that 

For a ID spin system with translation and an 
on-site projective symmetry u{g), the symmetric 
ground state cannot be short-range correlated, if 
the projective symmetry u{g) correspond to a non- 
trivial elements in H 2 (G,C). 



by a TI MPS. Its fixed-point MPS also has an on-site 
projective unitary symmetry u(g). For a proper choice of 
block size n, we can make u(g) and u(g) to be the same 
type of projective representation described by uj sym G 
H 2 (G, C). For TI fixed-point MPS, we have w [fc ] = U[k-i] 
since M[ fe ] (g) = M[k-i] (s)(cf. appendix D). Thus uj sym = 
0, that is, the trivial element in H 2 (G, C). So, if uj sym ^ 
0, the ID TI state cannot be short-range correlated. In 
other words 

ID spin systems with translation and an on-site 
projective symmetry are always gapless or have de- 
generate ground states that break the symmetries. 

If the ground state of the ID spin system does not break 
the on-site symmetry and the translation symmetry, then 
ground state is not short-range correlated and is gapless. 
If the ground state of the ID spin system breaks the 
on-site symmetry or the translation symmetry, then the 
ground state is degenerate. 

As an application of the above result, we find that 

ID half-integer-spin systems with translation and 
the 5*0(3) spin rotation symmetry are always gap- 
less or have degenerate ground states. 

which agrees with the well known result of Ref. 22. 

To have a gapped TI ID state with an on-site symme- 
try, the symmetry must act linearly (ie not projectively). 
In this case, we can show that the total phase factor of 
the state ai,(g) breaks up into L ID representations a(g) 
(see appendix F) 

For ID spin systems of L sites with translation 
and an on-site symmetry G, a gapped state that do 
not break the two symmetries must transform as 

u^{g)®...®u^{g)\4>L) = [a{g)] L \H) (21) 
for all values of L that is large enough. 

Here vS°\g) is the linear representation of G acting 
on the physical states in each site and a(g) is an one- 
dimensional linear representation of G. 

Let us apply the above result to a boson system with 
p/q bosons per site. Here the bosons number is conserved 
and there is an U(l) symmetry. Certainly, the system is 
well defined only when the number of sites L has a form 
L = Jq (assuming p and q have no common factors). 
For such L, we find that a^g) = ct (g) J = a (g) L ' q , 
where oto(g) is the generating ID representation of the 
U(l) symmetry group. So eqn. (21) is not satisfied for 
some large L. Therefore, 

a ID state of conserved bosons with fractional 
bosons per site must be gapless, if the state does 
not break the U(l) and the translation symmetry. 



The reason is as follows. If a ID state with translation 
symmetry is short-range correlated, it can be represented 



In higher dimensions, the situation is very different. A 
2D state of conserved bosons with fractional bosons per 
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site can be gapped, and, at same time, does not break the 
U(l) and the translation symmetry. 2D fractional quan- 
tum Hall states of bosons on lattice provide examples for 
such kind of states. 

Also, by repeating the discussions in section V A 1 for 
{nil-block TI LU transformations on the ID TI MPS, 
we can show that 

For ID spin systems with only translation and 
an on-site linear symmetry G, all the phases of 
gapped states that do not break the two symmetries 
are classified by a pair (to, a) where ui G H 2 (G,C) 
label different types of projective representations of 
G and a label different ID representations of G. 

Here a(g) is an ID representation of G that appear in 
eqn. (21). The symmetric LU transformations cannot 
change ID representation a(g). So the different phases 
are also distinguished by the ID representations a of G. 

Here are a few concrete examples: If we choose the 
symmetry group to be G — Z n , we find 

For ID spin systems with only translation and on- 
site J- n symmetry, there are n phases for gapped 
states that do not break the two symmetries. 

This is because Z„ has no projective representations and 
has n different ID representations. As an example, con- 
sider the following model 

# = D-^-^-i*K+i] ) ( 22 ) 

i 

where a x ' v ' z are the Pauli matrices. The model has a 
Z 2 symmetry generated by o~ z . The two different Z 2 
symmetric phases correspond the h —> 00 phase and the 
h —> —00 phase of the model. 

If we choose the symmetry group to be G — 5*0(3), we 
find 

For ID integer-spin systems with only translation 
and 50(3) spin rotation symmetry, there are two 
phases for gapped states that do not break the two 
symmetries. 

This is because 50(3) has only one ID representation 
and £f 2 (50(3), C) = Z 2 . Such a result agrees with the 
well known result that the AKLT state 40 of spin-1 chain 
and the direct product state with spin-0 on each site rep- 
resent two different 50(3) symmetric TI phases. The 
AKLT state (and the related Haldane phase 21 ) has gap- 
less boundary spin-1/2 states 41-43 and non-trivial string 
orders, 44,45 which indicate that the AKLT state is really 
different from the spin-0 product state. Actually, the full 
symmetry of 50(3) can be relaxed to only the dihedral 
group D 2 (J. 2 x Z 2 ) of rotation by ir around x, y and z 
axis. As explained in appendix C, D2 has one non-trivial 
projective representation, to which the AKLT state corre- 
sponds. AKLT is different from the spin-0 product state 
as long as on-site D 2 symmetry is preserved. This is 
consistent with the result in Ref. 20,45. 



5. ID systems with translation and parity symmetries 

In this section, we will consider the case of parity sym- 
metry for translational invariant system. We define the 
parity operation P for a spin chain to be in general com- 
posed of two parts: Pi, exchange of sites n and — n\ P 2 , 
on-site unitary operation u^°> where (u^) 2 = J. 55 Sim- 
ilar to the previous discussion, P gets redefined as we 
renormalize the state until at fixed point P\ becomes the 
exchange of renormalized sites and P 2 becomes u(°°) on 
every site, (it^ 00 )) 2 = I. The fixed point matrices hence 
satisfies (the 00 label is dropped): 

lZ n > ■ ,ir A ',, =±M~ 1 A i , i rM (23) 

3 l j r 

for some invertible matrix M, where we have used that 
the ID representation of parity is either (1, 1) or (1, — 1). 
We label the two ID representations with a(P) = ±1. 
Here M satisfies M~ 1 M T = e lB . But M = (M T ) T = 
e 2%s M, therefore, e %e = ±1 and correspondingly M is 
either symmetric M — M T or antisymmetric M = —M T . 
We will label this sign factor as /3(P) = ±1. 

Solving this equation gives that u — a(P)v(u l (g> u r ), 
where v is the exchange operation of two virtual spins 
i l and V and u l ,u r act on i l ,i r respectively. (u l ) T = 
f3(P)u l and [u r ) T — (3(P)u r . It can then be shown that 
each entangled pair \EPk,k+i) must be symmetric un- 
der parity operations and satisfies u k ® u l k+1 \E P k+ i_ k ) — 
a(P)\EPk.k+i) ■ There are hence four different symmet- 
ric phases corresponding to ct{P) — ±1 and /3(P) = ±1. 
By enlarging the local Hilbert space, we can show sim- 
ilarly as before that fixed points within each class can 
be mapped from one to the other with the {n.i}-block TI 
LU transformation preserving the parity symmetry. On 
the other hand, fixed points in different classes can not 
be connected without breaking the symmetries. There- 
fore, under the {ni}-block TI LU transformation, there 
are four block classes with parity symmetry and hence 
four parity symmetric TI phases: 

For ID spin systems with only translation and 
parity symmetry, there are four phases for gapped 
states that do not break the two symmetries. 

As an example, consider the following model 

H = J2l-BS! + Si-S i+1 }, (24) 

i 

where Si are the spin-1 operators. The model has a par- 
ity symmetry. The B — phase and the B — > +00 phase 
of the model correspond to two of the four phases dis- 
cussed above. The B = state 21 is in the same phase 
as the AKLT state. In the fixed-point state for such 
a phase, \EP kyk+1 ) = | fl) — | J/f). The parity trans- 
formation exchange the first and the second spin, and 
induces a minus sign: P : \EP k k+1 ) — > —\EP k k+ i). 
The B —> +00 state is the S z — 1 state. Its entangled 
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i±i a i±i a a i±i (a) 
i^u a a i^u a i±i (b) 
i^i a i^i a 1^1 a i^i (c) 
i^i a i^i a i^i a i^i (d) 



FIG. 4: (Color online) Representative states of the four par- 
ity symmetric phases, each corresponding to (a) a(P) = 1, 
P(P) = 1 (b) a(P) = -1, /3(P) = 1 (c) a(P) = -1, 
/3(P) = -1 (d) a(P) = 1, /3(P) = -1. + stands for a parity 
even entangled pair (e.g. 1 00) + |11}), — stands for a parity 
odd entangled pair (e.g. |01) — |10)). Each site contains four 
virtual spins. 



pairs are \EPk,k+i) = I tt) which do not change sign 
under the parity transformation. Thus the stability of 
the Haldane/AKLT state is also protected by the parity 
symmetry. 18-20 

To understand why there are four parity symmetric 
phases instead of two (parity even/parity odd), we give 
four representative states in Fig. 4, one for each phase. 
Connected pair of black dots denotes an entangled pair. 
+ stands for a parity even pair, for example 1 00) + |11), 
and — stands for a parity odd pair, for example 1 01) — 1 10) . 
Each rectangle corresponds to one site, with four virtual 
spins on each site. The four states are all translational 
invariant. If the parity operation is defined to be ex- 
change of sites together with exchange of virtual spins 1 
and 4, 2 and 3 on each site, then states (a) and (d) are 
parity even while (b) and (c) are parity odd. But (a) and 
(d) (or (b) and (c)) are different parity even (odd) states 
and cannot be mapped to each other through local uni- 
tary transformations without breaking parity symmetry. 
Written in the matrix product representation, the matri- 
ces of the four states will transform with a(P) — ±1 and 
(3(P) — ±1 respectively. Therefore, the parity even/odd 
phase breaks into two smaller phases and there are in all 
four phases for parity symmetric systems. 



C. Time Reversal Symmetry 

Time reversal, unlike other symmetries, is represented 
by antiunitary operator T, which is equivalent to the 
complex conjugate operator K followed by a unitary op- 
erator U . T has two projective representations: one on 
integer spins with T 2 = I and the other on half-integer 
spins with T 2 = —I. The classification of gapped ID 
time reversal invariant phases follows closely the cases 
discusses before. In this section, we will highlight the 
differences and give our conclusion. 



First, a state \<j>) is called time reversal invariant if 
T®T...®T\<j>) = fi\4>), where |/3| = 1. But for anti-unitary 
T, the global phase (3 is arbitrary and in particular we can 
redefine \<j>') = vW), such that T ® T... ® T\(jjf) = \<f>'). 
Therefore, in the following discussion, we will assume 
WLOG that (3 = 1. 

Now let us consider system without translational in- 
variance. T 2 = I or — / does not make a difference here 
as we can take block size 2 so that on the renormalized 
site, T 2 is always equal to /. Using argument similar to 
the case of on-site unitary symmetry, we can keep track 
and redefine symmetry operations as we do renormaliza- 
tion. Finally, at the fixed point we have a state described 
by matrices (j4^)(°°) which is invariant under time re- 
versal operation (T^ 00 ) = {u^)^K, that is, 



/ j i l i r .j l j 

j l i r 



N[k] = M[k-i] 



(25) 



where the fixed-point label oo has been omitted. Solving 
this equation we find, 

(a) M [fc] M [ * fe] = e ie I. As M [k] is invertible, e w = ±1. 

(b) w[ fc J = vW> 1 ® u^ r . where u^' l (u^ 1 )* = ±1 and 
u [k],r^ u [k],ry _ Therefore, each entangled pair is 
time reversal invariant 



(a 



[fc],r 



® U 



\ k + l U)K\EPi 



k,k+l) 



\EP, 



k,k+l/ 



(26) 



Similar to previous sections, we can show that uu* 
I and uu* = —I correspond to two equivalence classes 
and two time reversal invariant fixed point states can be 
mapped into each other if and only if they belong to the 
same class. Therefore, our classification result for time 
reversal symmetry is 



For ID gapped spin systems with ONLY time re- 
versal symmetry, there are two phases that do not 
break the symmetry. 



If the system has additional translation symmetry, we 
can similarly classify the TI phases and find that 



For ID systems with only translation and time 
reversal symmetry T, there are two gapped phases 
that do not break the two symmetries, if on 
each site the time reversal transformation satisfies 
T 2 = I. 



ID integer-spin systems are examples of this case. The 
Haldane/AKLT state and the spin product state are 
representatives of the two phases respectively 15,19,20 . We 
also have 



ID systems with translation and time reversal 
symmetry are always gapless or have degenerate 
ground states, if on each site the time reversal 
transformation satisfies T 2 = — I . 



ID half-integer-spin systems are examples of such case. 
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FIG. 5: Tensor-network - a graphic representation of the 
tensor-product wave function (27) on a 2D square lattice. The 
indices on the links are summed over. 
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VI. GENERALIZATION TO HIGHER 
DIMENSIONS 

In the last few sections, we classified symmetry pro- 
tected topological orders in one dimension, using (sym- 
metric) LU transformations. Can we use (symmetric) LU 
transformations to classify (symmetry protected) topo- 
logical orders in higher dimensions? 

In higher dimensions, the situation is much more com- 
plicated. First, infinitely many kinds of non-trivial topo- 
logical orders exist for class of systems without any 
symmetries. 7,48 A partial classification is given in Ref. 8 
for such a case in 2D. In the presence of symmetry, the 
phase diagram is even more complicated. Ref. 49 stud- 
ied 2D topological orders protected by time reversal and 
point group symmetry. So far, we do not have a detailed 
understanding of topological phases in the presence of 
symmetry. 

However, using similar arguments as that used for ID 
systems, we can obtain some simple partial results for 
higher dimensions. For example, we have 



For d- dimensional spin systems with only transla- 
tion and an on-site symmetry G which is realized 
linearly, the object (a, u)x, 002, cjd) label distinct 
gapped quantum phases that do not break the two 
symmetries. Here a labels the different ID repre- 
sentations of G and uji £ H 2 (G,C) label the dif- 
ferent types of projective representations ofG. 



Let us illustrate the above result in 2D by considering a 
tensor network state (TNS) on a square lattice, where the 
physical states living on each vertex i are labeled by rrij. 
A translation invariant TNS is defined by the following 
expression for the many-body wave function f ({mj): 



ijkl--- 



Ami Am 2 Am 3 ATn 4 
^ejif^jhkg^lqrk^tlsi ' 



(27) 



Here A™^ kl is a complex tensor with one physical index 
rtii and four inner indices i,j,k,l. The physical index 
runs over the number of physical states d on each vertex, 
and the inner indices runs over D values. The TNS can 
be represented graphically as in Fig. 5. 



FIG. 6: (Color online) (cji = 0, 1;oj2 = 0, 1) label four dis- 
tinct states in integer spin systems with translation and spin 
rotation symmetries: (a) (wi,W2) = (0,0), (b) (wi,^) = 
(0,1), (c) (wi,w 2 ) = (1,0), (d) (w l5 W2) = (1,1). The open 
dots represent 0-spins. The solid dots represent 1/2-spins. 
Two solid dots connected by a line represent a spin-singlet 
formed by two 1/2-spins. The state in (a) has an spin-0 spin 
on each site. The state in (b,c) has two spin-1/2 spins on 
each site (or 4 states per site that form spin 0+1 representa- 
tion of SO(3)). The state in (d) has four spin-1/2 spins on 
each site (or 16 states per site that form spin 0+0+1+1+1+2 
representation of SO(3)). 



If the tensor A satisfy 



A% ud = a(g) E u mm ,(g)x (28) 

I'r'u'd' 

{M~ l {g))u'M r/r (g)(N- 1 (g)) (id/ N ulu {g)Af' rlu , d ,, 



for some invertible matrices M(g) and N(g), then the 
many-body wave function ^({rrij}) is symmetric under 
the on-site symmetry transformation g in the on-site 
symmetry group G. Here a(g) is a one dimensional rep- 
resentation of G, the D by D matrices M(g) form a pro- 
jective representation represented by u>i £ H 2 (G, C), and 
the D by D matrices N(g) form a projective representa- 
tion represented by ui 2 £ H 2 (G,C). Since symmetric LU 
transformation cannot change (a, u>\, U2), thus (a, oji, uj 2 ) 
label distinct quantum phases. 

In fact (a, ui%, UJ2) are all measurable, so they indeed 
label distinct quantum phases. On a torus of size L x x 
L y , the symmetric many-body wave function $({rai}) 
transforms as the ID representation a LxLy (g) under the 
on-site symmetry transformation g. If G has only a finite 
number of ID representations, we can always choose L x 
and L y such that a L " Ly (g) = a{g). 

On a cylinder of size L x x L y with open bound- 
ary in the a;-direction, the states on one boundary will 
form a projective representation which is represented 
by L y uji £ H 2 (G,C). Similarly, if the open bound- 
ary in the y-direction, the states on one boundary will 
form a projective representation which is represented by 
L x u> 2 £ H 2 (G,C). If we choose L x and L y properly (for 
example to make L x u>2 = uj 2 and LyUi = wi), we can 
detect both u>i and uj2- 
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Symmetry 


No. or Label of Different Phases 


Example System 


None 


1 




On-site Linear Symmetry of Group G 


uj e H 2 (G,C) (*) 


On-site Z n or SU(2): 1 phase 
On-site SO(3) or D2 on integer spin: 2 phases 


On-site Projective Symmetry of Group G 


lj e H 2 (G,C) 


On-site 50(3) or D 2 on 
half-integer spin: 2 phases 


Time Reversal(TR) 


2 




Translational Invariance(TI) 


1 




Tl+On-sitc 
Linear Symmetry of Group G 


uj e H 2 {G,C) and a(G) 


TI+On-site Z n : n phases 
Tl+On-site SO(3) on integer spin: 2 phases 


TI+ On-site 
Projective Symmetry of Group G 





TI+On-site bO(3) or £>2 on 
half-integer spin: no gapped phase 


TI+Parity 


4 




TI+TR 


2 if T ' 1 = I 
if T 2 = -I 


TI+TR on integer spin: 2 phases 
TI+TR on half-integer spin: no gapped phase 



TABLE I: Summary of classification result for ID gapped spin system with symmetric ground states. TI stands for translational 
invariance. TR stands for time reversal symmetry. H 2 (G, C) is the second cohomology group of group G over complex number 
C. a(G) is a ID representation of G. (*): this result applies when a(G) form a finite group, when G = U(l), further 
classification according to different a(U(l)) exist. 



A system with integer on-site spins gives us an ex- 
ample with G = 50(3), if the system has translation 
and spin rotation symmetry. For G = 50(3) there is no 
non-trivial ID representations. So we can drop the con- 
sideration of a. Also G = 50(3) has two types of pro- 
jective representations: H 2 (SO(3),C) = {0,1}, where 
uj = is the trivial projective representation which cor- 
respond to linear representations of 50(3) (on integer 
spins) and co = 1 is the non-trivial projective repre- 
sentation which corresponds to half-integer spins. Thus 
{uj\ = 0, 1; uj 2 = 0, 1) label four distinct states in 2D (see 
Fig. 6). 



Results similar to ours on the classification of 
integer/half-integer ID spin chains with SO (3) symme- 
try have been obtained independently by Alexei Kitaev 
recently. 50 

We would like to thank A. Kitaev, J. I. Cirac, F. Ver- 
straete, G. Vidal, G.-M. Zhang, A. W.W. Ludwig and 
H. Katsura for some very helpful discussions. XGW is 
supported by NSF Grant No. DMR-1005541. ZCG is 
supported in part by the NSF Grant No. NSFPHY05- 
51164. 



VII. CONCLUSION 

Using the (symmetric) Local Unitary equivalence rela- 
tion between gapped ground states in the same phase and 
the matrix product representation of ID states, we clas- 
sify possible quantum phases for strongly interacting ID 
spin system with certain symmetry when ground state of 
the system does not break the symmetry. Our result is 
summarized in table I. 

Many well known results are rederived using a quite 
different approach, for example the existence of Hal- 
dane phase 21 and the gaplessness of spin 1/2 Heisenberg 
model 22 . Those results are also greatly generalized to 
new situations, for example to the case of time reversal 
symmetry and D 2 symmetry. We find that the projective 
representations play a very important role in understand- 
ing and formulating those generalized results. 

In higher dimensions, things are more complicated. 
Nevertheless, similar considerations allow us to obtain 
some interesting examples of symmetry protected topo- 
logical orders. A complete classification of higher dimen- 
sion phases, however, requires at least a full understand- 
ing of topological orders, an element that is absent in the 
ID phase diagram. 



Appendix A: Local Unitary Transformation on 
Matrix Product State and Invariance of Double 
Tensor 

In this section we present the proof for the theorem 
that: Two sets of matrices Ai iCc p and Bj^ a p, where i,j 
label different matrices and a/3 are the row and column 
indices of the matrices, give rise to the same double ten- 
sor §L al ,l3 X = T,i A t-aP X A *i, lx = J2j B 3,aP X B *j,-y X > ^ 

and only if they are related by a unitary transformation 
Bi, a p = ^2jUijAj^ a fj. The dimension M of i and the 
dimension N of j can be different. The M x N matrix 
U is in general called unitary if WU = UW = II, the 
identity matrix in L dimensional space where L < M and 
L < N. The complete proof of this theorem can be found 
in Ref. 35, where this property is discussed in terms of 
'the unitary degree of freedom in the operator sum rep- 
resentation of quantum channels'. Here we rc-prescnt 
this proof following the notation and terminology of the 
current paper for simplicity of understanding. 

First we prove the 'if part of the theorem. Suppose 
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that Bi^p = J2J u i3 A j,afs, then 

i 

i ji 32 

jl 

Therefore the two double tensors are the same. This 
proves the first part of the theorem. 

On the other hand, suppose that the two double ten- 
sors are the same E^ 7(3x = Ef 7( g x = E Q7S( g x . Reorder 
the indices of E and treat it as a matrix with row indices 
aj3 and column indices 7%. We will denote the double 
tensor after this reordering as Eq,^ j7X . It is easy to see 
that E is a positive semi-definite matrix, as 

X v * a ^nx v ix ( A2 ) 

i a/3 7X 

for any vector v a p 
Diagonalize E into 

E a ^, 7X = ^Afcefe !a) g x e£ )7X (A3) 

k 

with Xk > 0. Define vectors ek, a /3 = \/Xke-k,ap, so that 
Ea/3,7x = J2k^k, a /3 x e^. e fc , Q/3 form a complete or- 
thogonal set and hence we can expand Ai >a /j and Bi iCt p 
in terms of them. 

A,a/3 = X Pik&k,af3 (A4) 
-Bi,Q0 = ^ Qifce/c,q/3 

Then 

E Q /3, 7X = ^ ' Ai^ a p x ^4i i7X (A5) 

i 

= X X B ik 1 Pik 2 ^k u al3 X e£ 2i7x 

But we know that E Q/ 3 j7X = ^ fc efe )C(| g x e* k . Therefore, 
J2i p ikiP*k 2 = 4ifc 2 and P is a unitary matrix. 

Similarly we can show that Q is a unitary matrix. 
Therefore A^ a p and Pj, Q/ 3 is related by a unitary trans- 
formation Uij where U — PQ^ . We have thus proved 
both part of the theorem. 



Appendix B: Degeneracy of Largest Eigenvalue of 
Double Tensor and Correlation Length of Matrix 
Product State 

In section III, we cited the property of matrix prod- 
uct state that the finite correlation length of the state is 
closely related to the non-degeneracy of the largest eigen- 
value of the double tensor, as discussed in Ref. 29,30. In 
this section, we give a brief illustration of why this is so. 
For simplicity of notation, we focus on the translational 
invariant case first. Generalization to matrix product 
states without translational invariance is straightforward 
and similar conclusions can be reached. 

For a matrix product state |</>) described by matrices 
Ai, a p with double tensor E ai ,p x = J2i A i,a/3 X A* 7X > 
define E[0] a7i/3x = J^y OijAi >a p x A* j ryx for arbitrary 
operator Oy. Follow the previous convention and treat 
E and E[0] as matrices with row index a"f and column 
index fy- The norm of the wavefunction is given by 
(4>\4>) = Tr(E N ), where N is total length of the chain. 
WLOG, we can set the largest eigenvalue of E to be 1 and 
hence the norm goes to a finite value (dimension of the 
eigenspace) as N goes to infinity. The expectation value 
of any local operator O is (O) = Tr(E Ar - 1 E[0])/Tr(E JV ) 
and the correlation between two operators 0\ and 2 
becomes 

(Oi0 2 > - (Oi)(0 2 ) = 

Tr(E Ar - i - 2 E[0 1 ]E L E[0 2 ])/Tr(E JV )- (Bl) 
Tr(E Ar - 1 E[0 1 ])Tr(E Ar - 1 E[0 2 ])/Tr 2 (E JV ) 

The physical constraints on the expectation value and 
correlation functions of local operators require that the 
double tensor E has certain properties. First, put E 
into its Jordan normal form and decompose it as E = 
J2\ AP\ + R\, where P\ is the diagonal part and R\ 
the Nilpotent part. But for the largest eigenvalue 1, 
Ri must be as otherwise for large system size N 
(O) = Tr(E Ar " 1 E[0])/Tr(E Ar ) ~ Tt((Pl + P i ) JV E[0]) will 
be unbounded for any E[0] that satisfies Tr(i?iE[0]) 7^ 0. 
The physical requirement that any local operator has 
bounded norm requires that R\ must be 0. 

Next we will show that the dimension of Pi 
is closely related to the correlation length of 
the state. At large system size N, the correla- 
tor (O1O2) - (Oi)(0 2 ) = Hr(i\E[Oi](X; A Ai^ + 
i? A ) i E[0 2 ])/Tr(P 1 ) - Tr(P 1 E[0 1 ])Tr(P 1 E[0 2 ])/Tr 2 (P 1 ). 
When L is large, we keep only the first or- 
der term in AP\ + R\) L and the cor- 
relator goes to Tr(P 1 E[0 1 ]P 1 E[0 2 ])/Tr(P 1 ) - 
Tr(P 1 E[0 1 ])Tr(P 1 E[0 2 ])/Tr 2 (P 1 ). If P 1 is 
one dimensional, the two terms both become 
(vi|E[Oi]|t>i)(«i|E[02]|vi) and cancel each other for any 
Oi,0 2 and the second order term in (J^ A AP\ + R\) L 
dominates which decays as A L . For A < 1, the cor- 
relator goes to zero exponentially and the matrix 
product state as finite correlation length. On the 
other hand, if P\ is more than one dimensional, the 
first order term has a finite contribution indepen- 
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dent of L Etj<«i|E[0 1 ]|«i>(«i|E[0 3 ]|« i )/r&(Pi) - 

(vi\E[0 1 ]\vi)(v j \E[0 2 ]\v j )/'It i (P l ) i where v it Vj are 
eigenbasis for P\. Therefore, degeneracy of the largest 
eigenvalue of the double tensor implies non-decaying 
correlation. To describe quantum states with finite 
correlation length, the double tensor must have a largest 
eigenvalue which is non-degenerate. 

If the system is not translational invariant and E^ 
vary from site to site, we cannot diagonalize all the 
double tensors at the same time. However, as shown 
in the 'canonical form' of Ref. 30, there is a largest 
eigenspace of E^ (with eigenvalue 1) such that the right 
eigenvector on site k is the same as the left eigenvector 
on site k + 1. Therefore, when multiplied together this 
eigenspace will always have eigenvalue 1. There could be 
other eigenspaces with eigenvalue 1 and matching eigen- 
vectors from site to site. However, then we can show 
similar to the TI case that this leads to an infinite corre- 
lation length. On the other hand, other eigenspaces could 
have eigenvalues smaller than 1 or they have mis-matched 
eigenvectors. If this is the case, all other eigenspaces de- 
cay exponentially with the number of sites multiplied to- 
gether which gives rise to a finite correlation length. We 
will say that E^l has a non-degenerate largest eigenvalue 
1 for this case in general. 



Appendix C: Projective Representation 



an Abelian group, which is called the second cohomology 
group of G and denoted as H 2 (G,C). The identity ele- 
ment Wo of the group is the class that contains the linear 
representation of the group. 

Here are some simple examples: 

(a) cyclic groups Z n do not have non-trivial projective 
representation. Hence for G — Z n , H 2 (G,C) contains 
only the identity element. 

(b) a simple group with non-trivial projective represen- 
tation is the Abelian dihedral group D 2 = Zi x Z 2 , For 
the four elements of the group (0/1,0/1), consider repre- 
sentation with Pauli matrices g(0, 0) 



7(1,0) 



1 

-1 



5(1,1) 



1 

-i 

1 



.9(0, 1) 



It can be 



check that this gives a non-trivial projective representa- 
tion of D 2 . 

(c) when G = SO(3), H 2 (G,C) = Z 2 . The two elements 
correspond to integer and half-integer representations of 
50(3) respectively. 

(d) when G = 17(1), H 2 (G,C) is trivial: H 2 (U(1),C) = 
Z\. We note that {e lm6 } form a representation of 
U(l) — {e l9 } when m is an integer. But {e lm6 } will 
form a projective representation of U(l) when m is not an 
integer. But under the equivalence relation (C3), {e lme } 
correspond to the trivial projective representation, if we 
choose /3(g) — e~ lm6 . Note that /3(g) can be a discon- 
tinuous function over the group manifold. 



Operators u(g) form a projective representation of 
symmetry group G if 

u(gi)u(g 2 ) = u(gi,g2)u(gig 2 ), 51,52 G G. (CI) 

Here uj(gi,g2)in£., the factor system of the projective 
representation, satisfies 

^(52,53)^(51,5253) = u(gi, 52M5152, 53), (C2) 

for all 51,52,53 G G. If U)(g\,g 2 ) = 1, this reduces to the 
usual linear representation of G. 

A different choice of pre-factor for the representation 
matrices u'(g) = (3(g)u(g) will lead to a different factor 
system ui'(g 1 ,g 2 ): 

W(5l ' 52) = «« {9U92) - (C3) 

We regard u'(g) and u(g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems 0/(51,52) and 01(51,52) as be- 
longing to the same class uj. 

Suppose that we have one projective representation 
1x1(5) with factor system uii(gi,g 2 ) of class uii and an- 
other 112(g) with factor system ^2(51,52) of class iv 2 , ob- 
viously u\ (g)®u 2 (g) is a projective presentation with fac- 
tor group uji (51 , 52 )u>2 (51 , 52 ) ■ The corresponding class uj 
can be written as a sum u>\ + cj 2 • Under such an addi- 
tion rule, the equivalence classes of factor systems form 



Appendix D: Solving symmetry condition for fixed 
point 

In this section, we explicitly solve the symmetry condi- 
tion eqn. (14). The goal is to 1. classify possible symme- 
try operations at fixed point and 2. find the correspond- 
ing symmetric fixed point state. For simplicity, we drop 
the site index [k] and rewrite eqn. (14) as 

£ U ZSf(g) ) A ^ = ^(9)^^(9) (Dl) 

j l j r 

where u(g) is a projective or linear unitary representation 
of G, the matrix is given by its matrix elements 

A iH r jCt p = <J*\,6il a ■ y/X£6irp with l\ot = l,...,D t , 

i r ,(3 = l,...,D r , and M(g), N(g) are sets of invertible 

matrices labeled by 5. Since a (% g ^ is also a projective or 

linear unitary representation of G, we can absorb a^ R \g) 
into u(g) and rewrite eqn. (Dl) as 

J^Uii^rfoJA^, =N- 1 (g)A iH rM(g) (D2) 

j l j r 

We note that matrix elements A^ ir a p is non-zero only 
when a — i l , (3 = i r and the full set of {A^r} form 
a complete basis in the space of Di x D r dimensional 
matrices. 
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M(g), N(g) do not necessarily form a representation 
of G. But the fixed point form of the matrices requires 
that M(g), N(g) be a projective representation, as on 
the one hand 



u i'i r ,j l j'-(9l92)Ajlj 



(D3) 



rr 



T3 r 



= ^2 UJ sym{9l,g2)u i l i r k l k r(gi)u k l k r j l : j r (g2)A 

jl jr k l k r 



^u: sym {g ll g 2 )N'"{g 2 )N~ L {g 1 )A lH ,M{g 1 )M{g 2 ) 

and on the other hand 

^ti^r^-r {g 1 g 2 )A j i 3 t (D4) 

-A r_1 (.9i52)^ l ^M(<?i52) 

Therefore 

w sym (5i,32)A r_1 (52)^ 1 (5i)^ l -^(.9i)M(ff2) 
= JV-Vfla^'irM^fli) (D5) 

for all i z i r . However, the set of matrices {Atf^} form 
a complete basis in the space of Di x £) r dimensional 
matrices. Therefore, 

Usym (91 , 92 ) N~ 1 (.92 ) JV~ 1 (ffl ) ® M( 5l ) M ( ff2 ) 

= N- 1 (g 1 g 2 )®M(g 1 g 2 ), (D6) 



and iV(<7) and M(g) form two projective representations 

(D7) 



N(gh)=u N (g,h)N(g)N(h), 
M(gh)=u> M (g,h)M(g)M(h), 



with |wjv(ffi.,52)| = |wAf(fli)p2)| = 1 and 

wm (51,32) 



Wsym(gi,g2) 



(D8) 



UN (9l, 9n) 

Let us rewrite eqn. (D2) as 

^(5)E«i»i^ J i>jr(ff)^i»iO^~ 1 (ff) = A'r (D9) 



T3 r 

We note that 



Nig^ifi'^j^i'Mi^rA^M-^g) = A, 

3 l j r 

where the matrices M and N are given by 



(D10) 



(Dill 



Since the set of matrices {Ajjjr} form a complete basis 
in the space of x D r dimensional matrices, we find 



Uii t r,i r (g) = (N i i(g)M i r ij r(g). 



(D12) 



Putting back the factor of a^ R \g), we find that 

umr^rig) = aW( g )(N-i). liil (g)M.r tj r(g). (D13) 



Appendix E: Equivalence Between Symmetric Fixed 
Point States 

From the solution in section D, we know that the fixed 
point state symmetric under linear on-site symmetry of 
group G takes the form 



>(°°) = \EP ia )\EP 2 ^)...\EP k>k+1 ) 



(El) 



where \EP kk+ i) is an entangled pair between the right 
virtual qubit on site k and the left virtual qubit on site 
k + l(see Fig. 2 upper layer). Each entangled pair is 
invariant under a linear symmetry transformation of the 
form u^' r (g)t»u^ k+ ^' l (g) 



.Ik], 



(g) <g> u^ l (g)\EP k ^ +1 ) = \EP Kk+1 ) (E2) 



But u^' r (g) or u^- k+1 ^' l (g) alone might not form a linear 
representation of G. They could in general be a projec- 
tive representation of G. If uM' r (g) is a projective rep- 
resentation corresponding to class u> in H 2 (G,C), then 
u [k+i\,l mus t correspond to class — ui. uj does not vary 
from site to site and labels a particular symmetric fixed 
point state. 

Now we will show that symmetric fixed point states 
with the same uj can be connected through symmetric 
LU transformations and hence belong to the same phase 
while those with different u) cannot and belong to differ- 
ent phases. 

First, suppose that two symmetric fixed point states 
10!) and \(j) 2 ) are related with the same u>, i.e. 

u[ k] ' r (g) <g> u [k+1] \g)\EP kM+1 ) 1 = \EP kM+1 ) 1 (E3) 
uf' r {g) ®u 2 k+1] '\g)\EP k , k+l ) 2 = \EP kM+1 ) 2 

where |-E-Pfc.fc+i)i(2) is an entangled pair of virtual spins 
on Hilbert space ® Tif^' 1 ■ u i(2) (flO ^ s a projec- 

tive representation of G corresponding to uj on and 
u ^x(2p' 1 (&) a projective representation corresponding to 
-uj on rt 1(2) . 

We can think of \EP ktk+1 )i and \EP kyk+ \) 2 as liv- 
ing together in a joint Hilbert space (H^ ' r © H 2 r ) <X> 
('H^ k+1 ^' 1 © H 2 k+1 ^' L ). The symmetry representation on 
this joint Hilbert space can be defined as 



u^~{g)®u^\g) 



(E4) 



-K Ur ( 9 ) © u^{g)) ® (u[ k+ ^(g) © 4 fc+ 
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FIG. 7: (Color online) Fixed point state related to projective 
representation of class uj before (upper) and after (lower) a 
local unitary operation on the shaded region that does not 
break the symmetry. White dots correspond to ojo, the iden- 
tity element in H 2 (G, C), black dots correspond to w and gray 
dots correspond to — u>. 

As u [ ^ r (g) and uf' r {g) (also u^ +1] \g) and u [ 2 k+1] \g)) 
both correspond to u) (—uj), their direct sum 
v}- k '' r (g)(w- k+1 '' l (g)) is also a projective representation 
corresponding to uj(—uj). Therefore, we have a linear 
representation of G on each site k, u^' l (g) (g) u^' r (g) 
and both \EP k ,k+i)i and \EP k ,k+i)2 are symmetric un- 
der U W- r (g) (E)u^ k+1 ^ l (g). 

Now we can perform a LU transformation on the 
joint Hilbcrt space and rotate continuously between 
\EP k , k +i)i and \EP kM+1 ) 2 . That is, 

U(6) = co S ( e -)I-isin(l)(\a}{b\ + \b){a\) (E5) 

where \a) = \EP k ,k+i)i, \b) = \EP k , k+1 ) 2 and 9 goes 
from to 7T. By doing this locally to each pair, we can 
map \4>i) to \(f>2) (and vice verse) with LU transforma- 
tions without breaking the on-site symmetry of group G. 
Therefore, \<fii) and |</>2) belong to the same phase if they 
are related with the same uj. 

On the other hand, if |0i) and 1 02) are related to w\ 
and CJ2 respectively, we will show that they cannot be 
connected by any LU transformation that does not break 
the symmetry. 

Suppose that uj\ is non-trivial, we start with \(f>i) 
and apply a local unitary operation U to a finite region 
(shaded in Fig. 7). \<jj\) is composed of invariant singlets 
of symmetry group G. If U does not break the symmetry, 
the resulting state should still be composed of singlets. 
The singlet pairs outside of the shaded region are not 
changed while those overlapping with the shaded region 
can take any possible structure after the operation U. 

No matter what the change is, the right virtual spin 
on the site to the left of the region corresponding to uj\ 
should form a singlet with some degrees of freedom in 
the region. As the singlet is invariant under a linear rep- 
resentation of G, these degrees of freedom in the region 
must form a projective representation of G corresponding 
to —ait. These degrees of freedom could live on one site 
or distribute over several sites. However, the sites only 
support linear representations. Therefore, there must be 
some remaining degrees of freedom on the same sites that 
correspond to oji. These remaining degrees of freedom 
must form singlets again with other degrees of freedoms 



in the region that correspond to —uj\. We can continue 
this argument until finally some degree of freedom in the 
region corresponding to — uj\ connect with the left virtual 
spin on the site to the right of the region corresponding 
to u>\ and form a singlet. 

In Fig. 7, we illustrate one possible structure of singlets 
after operation U. White dots correspond to w , the 
identity element in H 2 (G,C), black dots correspond to 
uj\ and gray dots correspond to —u)\. 

Therefore, we can see that no matter what the sym- 
metric LU operation might be on \4>\), singlet entangled 
pairs related to uj\ must connect head to tail and cover 
the whole length of the chain. In other words, we can 
not shrink a chain of singlet entangled pairs related to 
non-trivial uj\ continuously to a point or change it to uj 2 
by acting on it locally and without breaking the symme- 
try. Hence fixed point states with different uj cannot be 
related to each other with symmetric LU transformation 
and hence belong to different classes. 

Appendix F: A proof of eqn. (21) 

A gapped TI state can be represented by a uniform 
MPS. After R steps of {n.;}-block RG transformation, 
we obtain a MPS described by matrices (A^, )^ which 
is given by eqn. (11). To describe a state that does 
not break the on-site linear symmetry, here (A^r)^) 
is invariant (up to a phase) under u^ R \g) on each site. 
Therefore 16 , 

J2^,fr(9) A 3 'r =a [K) {g)M- 1 (g)A %H rM(g) (Fl) 

j'j r 

must be satisfied with some invertible matrix M{g). Here 
we have dropped the RG step label R (except in a™ (5)). 
Each coarse grained site is a combination of Y\i=i n i orig- 
inal lattice sites and (g) form a ID representation of 
G. 

So if the number of sites has a form L = QHiLi n ^ 
then qj,(j) in eqn. (13) will have a form 

a L (g) = [aW(g)f (F2) 

for any value of Q. Now let us choose Q — Yif=i n 'i where 
Hill n i ancl IliLi n i have no common factors. The to- 
tal system size becomes L = Y[f=i n iY[f=i n i- We can 
perform, instead, a R' step of {n^}-block RG transfor- 
mation, which leads to a ID representation a^ R \g). We 
find that a^{g) in eqn. (13) will have a form 

a L (g) = [aM (g)] Q ' (F3) 

where Q' = L/ < = Ilti «<• T ^us 

a L (g) = [a^(g)]^< = [a™ (g)]^ (F4) 
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Since YiiLi n i an d IliLi n 'i nave n0 common factors, there 
must exist a ID representation a(g) of G, such that 

a {R) (g) = K<?)] nf -"% ^ R '\g) = [a(g)]^<_ 

(F5) 

Now eqn. (F2) becomes 

^(ff) = Ks)] gnf = ini = Ks)] 1 (F6) 

which gives us eqn. (21). 

Appendix G: Equivalence between gapped TI 
b-phases and gapped TI phases 

We have been using the quantum circuit formulation 
to study TI systems and classify b-phases. As the quan- 
tum circuit explicitly breaks translational symmetry, it 
is possible that each b-phase contains several different TI 
phases. (On the other hand, states in different b-phases 
must belong to different TI phases.) In this section, we 
will show that each b-phase actually corresponds to a 
single TI phase by establishing a TI LU transformation 
between states in the same b-phase. We will use the time 
evolution formulation of LU transformation (2) and find 
a smooth path of gapped TI Hamiltonian whose adiabatic 
evolution connects two states within the same b-phasc. 

First, as an example, we consider the case of TI only 
and show that there is only one gapped TI phase. Each 
TI SRC MPS is described(up to local change of basis) 
by a double tensor E which has a non-degenerate largest 
eigenvalue 1. E can be written as 

E Q7 ,^ X = E° 7( g x + E^ 7) g x = A a7 A/3 X + E^ 7) g x (Gl) 

where A is the eigenvector of eigenvalue 1 and E' is of 
eigenvalue less than 1. In the 'canonical form' 30 , A Q7 = 
\ a 5 ai , X a > 0. Obviously, E° is a valid double tensor and 
represents a state in fixed point form. We can smoothly 
change E to E° by turning down the E' term to from 
t = to t = T as 

E(t) = E° + (l- ~)E' (G2) 

Every E(t) represents a TI SRC MPS state. To see this, 
note that if we recombine the indices a/3 as row index 
and 7x as column index and denote the new matrix as 
E, then both E and E° are positive semidefinite matrices. 
But then every E(t) is also positive semidefinite, as for 
any vector \v) 

(v\E{t)\v) = (v\E°\v) + (1 - ±)(v\E>\v) 

= (l-±)(v\E\v) + ±{v\£?\v) >0 

E(t) is hence a valid double tensor and the state repre- 
sented can be determined by decomposing E(t) back into 
matrices Ai(t). Such a decomposition is not unique but 
WLOG, we can fix the decomposition scheme, so that 



Ai(t) vary continuously with time and reach the fixed 
point form at t = T(up to local change of basis). The 
state represented | </>(£)) hence also changes smoothly with 
t and has a finite correlation length as all eigenvalues of 
E(t) expect for 1 are diminishing with t. Therefore, E(i) 
represents a smooth path in TI SRC MPS that connects 
any state to a fixed point state (up to local change of ba- 
sis). 

How do we know that no phase transition happens 
along the path? This is because for every state \4>(t)), we 
can find a parent Hamiltonian which changes smoothly 
with t and has the state as a unique gapped ground 
state. 47 Following the construction in Ref. 29,30, we 
choose a sufficiently large but finite I and set the par- 
ent Hamiltonian to be H{t) = — h(t)k,k+h where 
h{t)k t k+i is the projection onto the support space of the 
reduced density matrix on site k to k + I at time t. Note 
that this Hamiltonian is translation invariant. For large 
enough I, h(t)k,k+l will always be D x D dimensional. As 
the state changes continuously, its reduced density matri- 
ces of site k to k+l changes smoothly. Because the dimen- 
sion of the space does not change, h(t)k,k+l also changes 
smoothly with time. Moreover, it can be shown that 
H{t) is always gapped as the second largest eigenvalue of 
E(t) never approaches 1 29 > 30 . Therefore, by evolving the 
Hamiltonian adiabatically from t — to t — T, we obtain 
a local unitary transformation connecting any state to 
the fixed point form, and in particular without breaking 
the translation symmetry. 

Because any TI fixed point state can be disentangled 
into product state in a TI way, we find that all TI ID 
gapped ground states are in the same phase, if no other 
symmetries are required. 

If the system is TI and has on-site symmetry, we need 
to maintain the on-site symmetry while doing the smooth 
deformation. A TI SRC MPS which is symmetric under 
on-site symmetry of group G is described by matrices 
which satisfy 

UiMAj = a(g)M- 1 (g)AM(g) (G3) 

j 

for some invertible projective representation M(g). The 
double tensor E hence satisfies 

a'P'j'x' 

(G4) 

where the group element label g has been omitted. Being 
the non-degenerate one-dimensional eigenspacc of E, E° 
must be invariant under the same transformation, and so 
does E'. Therefore we have 

E %,Px= E M ««> M W'( M *)^>M* x ,El w 
Ea 7lto = E M «*' M W'( M %y M xx' E W,0'x> 
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Now we smoothly change the double tensor as in Eqn.G2. 
Evidently, the symmetry condition Eqn.G4 is satisfied for 
all t. 

Decompose E(t) back to matrices Ai(t) so that the rep- 
resented state \4>{t)} changes smoothly with time. De- 
note the symmetry transformed double tensor as EM( g )- 
As E M ( g )(i) = E(t), there must exist a unitary operator 
u(g)(t), such that 

ZiMWMt) = M-\g)A i (t)M(g) (G5) 

3 

where u(g)(t) is a linear representation of G. Redefine 
u(g)(t) = u(g)(t) x a(g), then 

52uij(g)(t)Mt) = a(g)M- 1 {g)Mt)M(g) (G6) 

3 

As Ai(t) is chosen to be continuous with time, from the 
above equation we can see that u(g)(t) is also continu- 
ous in time. On the other hand, u(g)(t) form a linear 
representation of G. For all the cases we are interested 
in, the linear representations of G are discrete. There- 
fore, as u(g)(t) evolves smoothly with time, it cannot 
change from one representation to another but only from 
one equivalent form to another which differ by a unitary 
conjugation. That is, u(g)(t) = V '(t)u(g)V* (i) , with a 
continuous V(t). We can incorporate V(t) into the ma- 
trices Ai(t) and define A\{t) = J2j V ij(t)Aj(t), so that 
Ai(t) is symmetric under u(g) for all t. In the follow- 
ing discussion, we will assume that such a redefinition 
is made and the symmetry operation of the system will 
always be u(g) ® ... ® u(g). Therefore, the continuous 
evolution of E(t) from t — to t — T corresponds to 
a continuous evolution of short range correlated states 
|<^(£)) which is always symmetric under the same on-site 
symmetry u(g), with the same phase factor (a(g)) L and 
related to the same projective representation ui. 



Such a smooth path in symmetric state space cor- 
responds to a smooth path in symmetric Hamiltonian 
space. Construct parent Hamiltonian as discussed previ- 
ously. Because the state is symmetric under on-site u(g), 
the support space on site k to k + 1 must then form a rep- 
resentation space for (®u(g)) 1 . Therefore, it is easy to see 
that the parent Hamiltonian, being a summation of pro- 
jections onto such spaces, is also symmetry under on-site 
u(g). Moreover, the Hamiltonian remains gapped and TI. 
In this way, we have found a smooth path of symmetric, 
in particular TI, Hamiltonian whose adiabatic evolution 
connects any symmetric state labeled by a(g) and U) to 
the corresponding fixed point state (up to a local change 
of basis) and hence establishing the symmetric TI LU 
equivalence between them. 

As we show in appendix E that fixed point states with 
the same a(g) and u> can be related by symmetric local 
unitary transformations to each other, we now complete 
the proof that for ID spin systems with only translation 
and an on-site linear symmetry G, all gapped phases that 
do not break the two symmetries are classified by a pair 
(uj,a) where uj e H 2 (G, C) label different types of projec- 
tive representations of G and a label different ID repre- 
sentations of G. 

Similarly, if the system has translation and parity sym- 
metry, we can establish the equivalence between states 
labeled by the same a(P) and (3(P) in a translational 
invariant way (see the discussion below (23)). The pro- 
cedure is totally analogous to the on-site symmetry case, 
with the only difference that the symmetry condition for 
the matrices and double tensors become 

52j UijAj = ±M~ 1 A i M 

We find that for ID spin systems with only translation 
and parity symmetry, there are four gapped phases that 
do not break the two symmetries. 
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